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Abstract 

We study models where the superpartners of the ordinary particles have con- 
tinuous spectra rather than being discrete states, which can occur when the 
super symmetric standard model is coupled to an approximately conformal sec- 
tor. We show that when superpartners that are well into the continuum are 
produced at a collider they tend to have long decay chains that step their way 
down through the continuum, emitting many fairly soft standard model particles 
along the way, with a roughly spherical energy distribution in the center of mass 
frame. 



1 Introduction 



In this paper we take the first steps in studying the collider phenomenology of super- 
symmetric models that have some approximate conformal symmetry, which is reflected 
in the fact that the superpartners of the standard model (SM) particles have a con- 
tinuous spectrum [1] rather than being discrete states. The superpartners can have 
continuous spectra when the supersymmetric standard model is coupled to a confor- 
mal field theory (CFT) sector [1] which gives rise to so-called "unparticle" behavior!] 
A mass gap is generated for the superpartners of (and also the excitations of) the SM 
particles if the conformal symmetry is softly broken in the infrared (IR). This scenario 
is most easily modeled in a five-dimensional (5D) anti-de Sitter (AdS) space using the 
AdS/CFT correspondence [6j[7]. The 5D AdS space is cut off by an ultraviolet (UV) 
brane where supersymmetry (SUSY) is broken. In the absence of an IR brane, the SM 
fields propagating in the bulk will have continuous spectra. Such a theory would have 
been ruled out if the spectra continue down to zero mass. However, a mass gap for 
nonzero-modes can be generated by introducing a soft wall in the IR, which can be 
parameterized by a position dependent bulk mass term or a dilaton field with an ap- 
propriate profile. The boundary conditions on the UV brane remove half of the fermion 
zero-modes so that the SM chiral fermions can be obtained. In the supersymmetric 
limit, the 4D theory consists of SM particles and their superpartners as the zero- modes, 
plus a continuum of Kaluza-Klein (KK) excitations starting from some gap for each 
field. After SUSY breaking is introduced on the UV brane, the zero modes of the su- 
perpartners are lifted while the mass gaps of the continuum excitations are not affected 
since they are determined in the IR. Depending on the parameters, for large enough 
SUSY breaking the zero mode of the superpartner can merge into the continuum and 
only a continuous spectrum for the superpartner is left. 

As one can imagine, the collider phenomenology could be quite complicated with 
continuous spectra. Calculations of production cross-sections have already been dis- 

1 Unparticles are fields with a continuous spectrum [HG2, possibly with a with a mass gap [H[S], 
whose two point functions exhibit a nontrivial scaling behavior. 
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cussed in [5] , so here we will focus on decay processes. Imagine that some highly excited 
mode in the gluino continuum is produced at the Large Hadron Collider (LHC), it can 
decay to a squark with an arbitrary mass in the continuum between the squark gap 
and the initial gluino mass (neglecting the mass of the emitted quark). The squark can 
then decay back to a gluino as long as the squark mass is above the gluino gap. An 
obvious question is: does a continuum decay prefer to occur for small mass differences 
or large mass differences? If the gluino prefers to decay to a squark with an invariant 
mass close to its own, then the jet emitted in the decay will carry a relatively small 
amount of energy, and there will be many steps of decays before it reaches the bottom 
of the spectrum. The events will contain a high multiplicity of soft visible particles, 
which can be quite challenging at hadron colliders. On the other hand, if the gluino 
prefers to decay to the squark near the bottom of the spectrum, then we expect only a 
few decay steps and hard jets from the decays. The signal events in this case are more 
like traditional SUSY models. Since the theory becomes conformal at high energies, we 
expect that in the high invariant mass limit we will be closer to the first picture [8HTT] . 
Here we would like to make some more quantitative statements. 

Explicit calculations for continuum superpartner events, however, are not straight- 
forward. Because the initial and/or final states are not particles, the usual Feynman 
rules for particles are not directly applicable. One way to avoid this problem is to 
introduce a regularizing IR brane which makes the extra dimension compact, then the 
continuum becomes discrete KK modes [12] and we can perform calculations just as 
in the particle case. The continuum limit is obtained by taking the position of the IR 
brane to infinity. The physical results should not depend on this position, as long as 
it is much larger than the length scale of the inverse mass gap so that the KK modes 
are dense enough to approximate a continuum. Even in this case, one may worry 
about whether the usual narrow-width approximation of splitting a decay chain into 
steps with independent decays is a good one, as there are an infinite number of KK 
modes which can make virtual ("off-shell") contributions, reflecting the fact that the 
continuum state do not have a mass shell to be on. We will discuss the validity of this 
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approach and compare full calculations with the narrow-width approximation. 

This paper is organized as follows. We first review the 5D construction of con- 
tinuum superpartners, then plough into the details of the decay chains, showing that 
the narrow-width approximation is generally valid in perturbative theories. We then 
discuss the phenomenology and the parametric dependence of the observable quantities 
and give conclusions. A detailed check of the IR regulator is included in the Appendix. 

2 A Review of Continuum Superpartners 

Let us recall the setup used in Ref. [I], which will be the starting point of our col- 
lider studies. We consider a 4D SUSY theory with approximate conformal symmetry. 
Through the AdS/CFT correspondence [6] this can be modeled by a 5D AdS space. 
We take the metric of the AdSs space to be 



The space is cutoff at z = zuv = e by a UV brane where SUSY is broken. As is well 
known, such a theory can be described in the language of 4D M = 2 superfields, which 
implies that for each matter field its 5D M = 1 hypermultiplet \I/ can be decomposed 
into two 4D M = 1 chiral superfields $ = {</>, x, F} and $ c = {0 C , tp,F c }, with the 
fermionic Weyl components forming a Dirac fermion |T3j . In the case of gauge fields, 
a 5D M = 1 vector superfield can be decomposed into an M = 1 4D vector superfield 
V = {A^, Ai, D} and a 4D M = 1 chiral superfield o = {(£ + iA 5 )/y/2, A 2 , F a }. 

As in the usual case of extra-dimensional theories, one can decompose the matter 
fields as: 



where the relationships between scalar and fermion profiles are provided by SUSY. 




(1) 
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With a fermion bulk mass m(z) = c + fifZ, the bulk equations of motion give 
° 2 * , ^2 ,2 o/^/c c(c- 1 



/s+ ^-^-2^-nr_^j /fl = o, (4) 
U+{/-fi-^- C ^)h = 0. (5) 



The solutions for 2) and z) can then be expressed as linear combinations of 
the first order and second order Whittaker functions: 



1 „ r~, : , , „„ 1 



f L (p,z) = a M(k,- + c,2 J fi-j?z) + bW(K ,- + c,2Jti*-p*z) , (6) 



2{l + 2c)J^-p 2 x 
fn{p,z) = -a j M(K,-- + c,2yJ/i 2 f -p 



■ 2 z) 



-b P W(k, -I + c, 2 Jim} - p 2 z) , (7) 

+ yV/ - V 2 ) 1 



Cfif 



flj-p 2 



where a and b are coefficients which may depend on four-momenta p and that are 
determined by boundary and normalization conditions. We will take all SM fields as 
left-handed, so we demand that the right-handed chiral superfield $ c (p, z) satisfies the 
Dirichlet boundary condition at z — Zuv-, Q c {p,zuv) — 0. The left-handed chiral 
superfield $, through the equations of motion, satisfies modified Neumann boundary 
conditions. As a result, only the left-handed chiral superfield has a normalizable zero- 
mode. As shown in Ref. PQ, the zero mode profile for the left-handed field is given 
by 

f° L (p,z)=Af(fi f ,0)z- c e-^, (9) 
where the normalization factor A/"(/x/, 0) is obtained from J °° /(0, z)f(0, z)*dz = 1, 

Af(n f , 0) = (2- 1+2 >7 1+2c r(i - 2c))- 1/2 . (10) 

In Fig. [TJ we show the zero mode profiles for two values of c. As can be seen, when c is 
positive, the zero mode is localized near the UV brane, while for c < 0, the zero mode 
is repelled away from the UV brane. 
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Figure 1: Normalized zero mode profiles. In the left panel we take c = —0.45, fif 
0.4 TeV while in the right panel, we take c = 0.45, and /i/ = 0.4 TeV. 



Similarly for gauge fields [T] we write: 



\i{p,z) 
A 2 (p,z) 



Zuv 



h L , 



zuv 



1/2 



h L , 



h 



R, 



Mpy 



z \ 



3/2 



■ 



h 



R 



111 



(12) 



y zuv J \ z uv, 
where h^R represents f^R evaluated at c = 1/2 and \x g is related to the dilaton vacuum 
expectation value [7], (<&) = e~ 2 ^ 9Z / g 5 . 

After including SUSY breaking on the UV brane, the zero mode of the superpartner 
will be lifted, and the SUSY breaking mass that it acquires depends on the overlap of 
its wave function with the UV brane. For c close to +1/2, the zero mode of the super- 
partner will acquire a SUSY breaking mass near the full strength (i.e., comparable to 
the SUSY breaking on the UV brane). On the other hand, if c < 0, the SUSY-breaking 
mass of the zero mode superpartner is suppressed relative to the SUSY breaking on 
the UV brane. Because c = 1/2 for the gauge field, we will take all SM fields to have 
c close to 1/2 in this paper in order to have similar SUSY-breaking masses. 

The spectrum of nonzero-modes is continuous without an IR cutoff. As mentioned 
in the Introduction, it is convenient to introduce a regulating IR brane at a large 
distance z = zjr = L so that we can deal with discrete normalizable KK states [T2] . 
The continuum limit is obtained by taking L — > +00. The coefficients a and b of the 
wave functions in Eq. (EUTJ) can be obtained by imposing the boundary condition on 
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the IR brane and the normalization condition, 

f L {p,z) = N L {Hf,p) \ M(k, ~ + c, 2yV/ - p 2 z) 
+ b . W ( K ,l + c ^ fjl *-p*z?) , 



(13) 



2(l + 2c)Jfi}-pi x 

f R (p,z) = Af R (fi f ,p)\ V .M(K,-- + c,2J$-p*z) (14) 

+ b-W(K,~ + c,2^}-p 2 z)\ . (15) 
Imposing the boundary condition fR(p,z uv ) = we determine that, 



M{k, -\ + c, 2Jn) - jfizuv) 2(1 + 2c) J y* - p 



(16) 



W(k, -§ + c, - p 2 ^c/y) Hf- ^jn)-p 2 

and /t = — c/i// y"j — p 2 . 

We are especially interested in the behavior of the superpartners in the conformal 
limit at high energies. In the limit of p S> /i/, the normalization factor can be expressed 
as 

The gluino profiles can be found similarly by making the replacements c — >■ 1/2 and 
jif — > /i 9 in Eqs. ( TTBTfTp]) . 

The KK spectrum is determined by the boundary condition at the IR brane, 
fn(p, zjr) = 0, which leads to: 

In the limits of interest, < p < l/^c/y, the left-hand-side of Eq.f lTH]) tends to 

zero, which implies that the numerator of the right-hand-side of Eq.f ll8p satisfies: 



M(-—0=^, 0, 2^1 -p* z IR ) oc cos(-7r - x /p* - y\z m ) = 0. (19) 
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Figure 2: Spectra for the gluino (blue dots) and the squark (red dots). The plot shows 
mass versus KK mode number. The parameters are chosen to be zyy = 10 -3 TeV -1 , 
zir = 80 TeV -1 , c = 0.5, fi g = 0.4 TeV and (Mf = 0.5 TeV for the left panel, and 
Zuv = 10- 3 TeV -1 , zir = 100 TeV" 1 , c = 0.5, fx g = 0.3 TeV and fx f = 0.4 TeV for the 
right panel. 

Solving Eq. ( fl9|) . we obtain an approximate expression for the KK masses, 

m l ~ V 2 g + {\ + n ) 2 v 2 /z 2 IR with n = 0, 1, 2, • • • . (20) 

Some sample spectra are shown in Fig. [21 where we have taken Zuy = R = 10 -3 TeV -1 
and c = 1/2. With mass gaps on the order of half a TeV, a 20 GeV mode spacing is 
quite a good approximation to a continuum. 

3 Continuum Decay Chains and Narrow- Width Ap- 
proximations 

Now we would like to study the decay of a continuum superpartner. The question 
is: what are the characteristic features of the decay process? Does it prefer to decay 
through multiple steps and emit soft particles during the decays, or to decay directly 
down to the bottom of the spectrum together with a hard particle? What is the typical 
energy of the visible particles produced in the decay chain relative to the energy and 
other parameters of the superpartner? As mentioned earlier, in order to use the calcu- 
lation techniques developed for particles, it is convenient to include a regularizing IR 



Figure 3: Feynman diagram for sequential KK gluino decay. 



brane to discretize the continuum into KK modes [12J, so that we can deal with ordi- 
nary particles in the initial and final states. To study the validity of the narrow-width 
approximation of the decay chain, which splits the process into a sequence of indepen- 
dent decays at each step, we compare the result from a calculation of a sequence of real 
two 2-body decays with that of the 3-body decay which includes all the contributions 
from the virtual intermediate continuum superpartner. We assume that we start with 
an initial state with an invariant mass much higher than the mass gaps and the SUSY 
breaking mass. To simplify the calculations, we will ignore the SUSY-breaking and the 
masses of the SM particles. We find that, in the case that the calculations are reliable 
(i.e., independent of zjr = L as long as L ^> p" 1 , fij 1 , fi^ 1 ) , the virtual contributions 
are indeed smaller than the real contributions so that the narrow-width approximation 
is reliable. 

3.1 The 2-body gluino decay calculation 

Once we include a regulating IR brane to make the spectra discrete, the calculation of a 
2-body decay should be straightforward. We consider an initial state of a gluino of KK 
level m, decaying in its rest frame to a quark and a KK squark of level n, g m — > u° L u*£ n ■ 
Since we ignore the quark mass, the decay can occur to any KK squark lighter than 
the initial gluino. The energy of the emitted quark depends on which level of the KK 
squark the gluino decays to. To calculate the decay rate, we need the coupling between 
the gluino, squark, and the quark. The effective interaction of the KK gluino, KK 
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squark and the quark in momentum space is given by: 

(2tt)" 4 [2txY U ° L(yPm ~ ^ ^(^n) 9 n (Pm) , (21) 

where c(p m , q n ) is the vertex coefficient which can be computed by integration of the 
gluino, squark and quark 5D profiles over the fifth dimension. The expression is: 



f ZlR / ZjJV\^ ( Z \ 

c(p m ,q n ) = N~ g (jj<g,p m ) Nu(jj, f ,q n y N u (jj,f,ti) g$ I dz [— J ( - — ) 



zuv 



2 

II r y 

z 



\ 3/2 / \ 2 

— h(qn,z)* e uz (—) h L (p m ,z), (22) 
zuv J \zuvj 

where g§ is the 5D gauge coupling, which only enters the calculation as an overall 

factor. The 4D gauge coupling is related to the 5D gauge coupling g 5 by integrating 

over the zero mode profiles in the kinetic term for the gauge fields coupled with the 

dilaton. In that case one finds j7] that, 

^2 = Jh 1 ^ (23) 

4 z uv log(l/ (2zuvfi g )) - 1e) ' 

where 7_e is Euler's constant. Since we are interested in the gluino decay into a quark 

and a squark, we consider the QCD gauge coupling which equals a s {m z ) = 0.1184 at 

the Z-mass. If we take zuv — 10~ 3 TeV -1 and the nominal value for the gluino mass 

the 5D strong gauge coupling assumes the dimensionful value 
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gap of Hg = 0.1 TeV, 
g 5 = 0.108 TeV~ 1/2 . 

An immediate observation from Eq. ( l2~2"j) is that a finite coupling in the limit zm — > 
+oo is obtained only for /// > fi g . For /i/ < fi g , the integral is dominated by the large 
z region and it blows up as we take zjr — > +oo. In this case we cannot perform a 
sensible calculation. One can see that the exponentially growing factor comes from 
the dilaton profile for the KK gluino, which is a result of that the zero mode gauge 
field has to be a constant along the fifth dimension due to gauge invariance. A similar 
situation happens in the Randall-Sundrum 2 (RS2) scenario [H]. If one calculates 
the self interactions among KK gravitons in RS2 with a regulating IR brane, one also 



2 The dependence on the gluino mass gap is only logarithmic, reflecting the gauge coupling's 
running. 



finds that the coupling blows up as one takes the IR brane to infinity. In the 4D 
CFT picture, these KK gravitons correspond to the conformal bound states. The large 
coupling just means that these bound states are strongly coupled. These divergent 
couplings can never make any physical process involving UV zero-modes infinite. As 
is well known, the KK picture sometimes can give misleading results when locality 
in the extra dimension is involved. The point is that any process starting with zero- 
modes localized on the UV brane will mostly be sensitive to the physics near the UV 
brane due to the locality in the extra dimension. The "nonlocal" process of producing 
KK gravitons in the deep IR region must be suppressed due to interference of various 
diagrams even though each of them can have a large coupling. In our theory, all zero- 
modes are localized near the UV brane (for c ~ 1/2), so all physical processes originated 
from the zero-modes should also only be sensitive to the physics near the UV brane. 
The divergent coupling in the /if < /i g case just means that the naive calculations 
done in the KK picture are not valid. On the other hand, for /if > jj, g the integral is 
dominated by the region near the UV brane, and the calculations are trustworthy. For 
the same reason, the coupling between a KK gluino, a KK squark, and a KK quark is 
also dominated by the IR region and diverges when the IR brane is moved to infinity. 
However, since the wave functions of the KK quarks are suppressed in the UV region, 
the decay to a KK squark plus a KK quark should also be suppressed if the initial KK 
gluino was produced in the UV region in the first place. To avoid such complications, 
we will restrict our study to the /if > fj, g case and consider decays to the quark zero 
mode only in this paper. 

After obtaining the coupling, it is straightforward to calculate the decay rate to 
each KK squark. To compare with the continuum limit, we express the result in terms 
of the differential decay rate as a function of the energy of the outgoing quark: 



dE u o ^ AE u o 

L nefi f <E n <E 



c{p m) qn)c\p m q n ) E \ Tr[t Q t a ] 

nfc iif<E n <Em L 



where t a are the SU(3) generators in the fundamental representation and p m = ^Jp 
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Figure 4: Gluino differential decay rate with respect to the first quark energy. The 
initial KK gluino mass is p g = 1.15 TeV and we use g§ = 0.108 TeV" 1/2 . We choose 
Zuv = 10- 3 TeV" 1 , z IR = 100 TeV -1 , c = 0.5, fi g = 0.3 TeV and /i f = 0.4 TeV. 

From conservation of energy-momentum we have, 



A typical differential decay rate as a function of the outgoing quark energy is shown 

in Fig. HI We have normalized the decay rate by a factor p ■ Zm as we do with all 

differential decay rate figures hereafter, in order to cancel the unphysical IR dependence 

Am n = tc/zjji in the continuum limit 1^1 At small quark energies, the decay rate is 

suppressed by phase space, while at large quark energies, the decay rate is suppressed 

by the couplings due to the small overlap between the gluino and squark wave functions 

of large KK level differences. Overall we see that the suppression due to the coupling 

3 It is easy to see that the coupling is proportional to L from the normalizations of the gluino 
and squark wave functions. The differential decay rate of a single KK gluino is proportional to 
{1/L) 2 x L = 1/L as the density of the final KK squarks is proportional to L. This is related to the 
fact that the overlap of a single KK gluino wave function with the UV region where the zero modes 
are located is also proportional to 1/L. In reality the initial KK gluino will be produced with some 
energy range and the number of KK gluinos in that energy range is again proportional to L, which 
cancels the unphysical 1/L dependence in the differential decay rate of a single KK mode. 




(25) 



which implies 




(26) 
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is stronger and the outgoing quark has a relatively soft spectrum. The suppression 
of high energy emissions is exactly the behavior that is expected in a CFT [8HTT]. 
Starting with very high KK modes this leads to approximately spherical events [TT] . 
A more detailed discussion of the quark spectrum and its parameter dependence will 
be presented in the next section after the narrow-width approximation is justified. 

3.2 The 5D mixed position-momentum propagator 

From the previous subsection we see that the result of the 2-body calculation shows 
that the gluino prefers to decay to the squark with a mass not far below the gluino mass. 
This means that the resulting squark is likely to decay again back to the gluino and 
the whole process will involve a long decay chain. An important question is whether 
the sequence of decays can be treated independently and correlations implied by the 
full intermediate propagator can be safely neglected. This question is less trivial for 
the continuum case than the usual particle case because there are an infinite number, a 
continuum, of intermediate states which can give contributions at each step. To study 
this problem, we consider the 3-body gluino decay process g m — > u° L u ^g n where the 
continuum squark is in the intermediate state. We will calculate the contribution from 
the virtual ( "off-shell" ) intermediate squark while using real KK modes for the initial 
and final states. 

In the KK picture with a regulating IR brane, the intermediate squark propagator 
is simply a sum of the particle propagators of all KK levels [12] . The propagator of an 
individual scalar particle is 

lim % - = 7i5(q 2 -m 2 n )+iV = 1 _ (27) 

e->o+ q z — m„ + le q z — m A n 

where V denotes the Cauchy principal value. The delta function (the real part) rep- 
resents the phase space for a real intermediate particle (the "on-shell" contribution), 
while the second term (the imaginary part) represents the contribution of a virtual 
intermediate state (the "off-shell" contribution). However, the KK picture is not the 
most convenient one to perform calculations with virtual intermediate states as it in- 
volves an infinite sum of singular functions. We know that in the continuum limit, 
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the infinite sum of KK propagators simply becomes the unparticle propagator [T2] . 
and the series of poles on the positive real axis of p 2 merge into a branch cut of the 
unparticle propagator. Therefore, we will employ the full unparticle propagator for 
the intermediate state in our calculation. In this subsection we derive the unparticle 
propagator in mixed position-momentum space [T5] . 

The left-handed scalar propagator satisfies the equation of motion for the left- 
handed profile, with a delta-function source, 

d 2 z + -d. + (c 2 + c - + + $ - p 2 ) l P{p, z, z>) = 3 5(z - 0,(28) 

z A z 2 z / \zuvj 

and the UV boundary condition for the propagator is 

1, 3 



+ Hfz) \ P(p,z,z') 



= 0. (29) 

z=z uv 

To connect to the continuum limit, with outgoing boundary conditions, we demand 
that the propagator is exponentially damped [15] at large Euclidean momenta: 

P(ip, z IR , z') \ P ^ +00)Z , <ZIR -> e- pzi *. (30) 

In order to solve Eq. (I28H291) for P(p, z, z'), we look for solutions in the regions z > z', 
P>(p, z, z'), and z < z', P<(p, z, z'), and use matching boundary conditions at z — z', 

P<(p,z,z , )-P > (p,z,z')\ lt=z , = 0, (31) 

\ 3 

Z 



d z P < (p,z,z!)-d g P > (p,z,J)\ g=s , = -i[ . (32) 

We can write the general solution to Eq. ( )28|) in terms of two independent solutions to 
the homogeneous equation [15], K{p,z) and S(p,z): 



X3/2 w(K,l + c,2Jfi 2 f -p 



- 2 z) 



K(p,z) = , (33) 

\zuvj W{k,\ + c, 2^ 2 - p 2 z uv ) 

S(p,z) = 



z \ 3/2 i r(i + c-«) 



z uvJ 2^ n}- j? r ( 2 + 2c ) 
l ~ l 



M(k, - + c, 2^/1/ " P 2 z) W(k, - + c, 2yV/ ~ p 2 z uv ) 
-W(k, l - + c, 2^ 2 -p 2 z) M(k, l - + c, 2^)-p 2 z uv Yj , (34) 
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which satisfy the following boundary conditions: 

K(p,z uv ) = l, K(ip,z)\ p _, +00 -»■ e~ v \ S(p,z uv ) = 0, S'(p, z uv ) = 1. (35) 

In the region z < z'(z > z'), the general solution can be written as 

P<(>)(p,z,z') = a <{>) K(p,z) + b <{>) S(p,z). (36) 

The boundary condition at the UV brane, Eq. (1291) . fixes the ratio of a < /b < to be 
proportional to the kinetic function Sf c (p) already encountered in Ref. [Tj: 



6< z uv P 2 

I 1 I — 

V /=P+^' 2 



(37) 



E Fc (p) = zuv- I — 1 -f ■ (38) 

The kinetic function Ep c (p) not only fixes the spectral density, but it is also essential in 
determining the phase space. Using the required IR behavior of P(p, z, z'), we conclude 
that P>{p, z, z') can only be proportional to K(p,z), (i.e., b > = 0). At this stage we 
use the matching conditions at z = z', so that in the range z < z', the left-handed 
squark propagator can be expressed as [131 : 



iP(p,z,z') = ^^-K{p,z)K{p,z')-S{p,z)K{p,z') . (39) 
zuv 

For the case z > z' we just interchange z •f-)- z' in Eq. ( 139]) . 



Obviously for p 2 > /x 2 , — p 2 becomes imaginary, so the propagator has a 

branch cut on the real axis for p 2 > fi 2 , and the discontinuity is just twice the real 
part. This discontinuity corresponds to a real intermediate unparticle (the "on-shell" 
contribution). From the analogy of the particle propagator (12 7p we interpret the imag- 
inary part of the unparticle propagator as the virtual ( "off-shell" ) contribution to the 
3-body decay process, while the real part (the discontinuity) corresponds to the phase 
space of the unparticle. In fact, one can use this phase space to calculate directly 
the real 2-body differential decay rate considered in the previous subsection, instead 
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of summing over KK modes in a discretized theory. In the Appendix we show the 
equivalence of the two pictures in the limit that the IR regulator is removed, L — > oo. 
The numerical results for the 2-body differential decay rates using the two approaches 
also agree well for large L. 

3.3 The 3-body gluino decay 

With the unparticle propagator derived in the previous subsection, we can compute 
the virtual contribution to the 3-body decay process, g m — > u L u°^g n . We take both the 
initial and final gluinos to be KK states in an IR regularized theory, but use the unpar- 
ticle propagator for the intermediate state with only the imaginary part corresponding 
to the virtual contribution. It is straightforward to find the amplitude squared for the 
virtual 3-body process by integrating the imaginary part of the propagator P(q, z, z') 
derived in Eq. (l39|) over the positions of the two vertices in the extra dimension, 

\M(g m -+ u° L u° L *~g n )\ 2 = Agt \v(p~ gm , q,p gn )\ 2 (p^.^o.) (p u o.^«) . (40) 

We have labeled the 4-momenta of the initial, final, and intermediate states as follows: 
p g ™ for the initial gluino g m , p g n for the outgoing gluino g n , p u o for the quark, p u o* for 
the anti-quark, and q for the intermediate squark. The factor v(p m ,q,p n ) is given by 
the integration in the extra dimension of the respective profiles and the propagator: 



v(plq, P f) = Af^ f ,0W g (fi g ,pfW g ^ g ,p^ g ) 

fZlR fZlR 

x dz dz' e^ fl f )z z l/2 ' c h L (p g n ,z)z^P(q,z,z') 

J Z UV J Zjjy 

e^~^ )z ' z ll l 2 - c h L (p^z') , (41) 

/ \ 3/2 / , n -3/2 

here we use a rescaled scalar propagator, P(q, z, z') = ( J ( \ P(q, z, z'), 
and we can write the differential decay rate as: 

4& = 9i v{p?, Q,P>(pT, Q, E ^ 2Eu ^ ~^ )2 \ ^) 2 ) 2 l Tr[t a t¥ta] (42) 
dE u y5 yiJ a^^g> g > 32pf(2E u -p q n )7r 3 8 1 



where _E M o is the energy of the (first) outgoing quark in the initial gluino rest frame 
and \Ti[t a t h t b t a ] = \C A C\ = 2/3 for SU{3). The range of £ u o is determined by the 
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Figure 5: Gluino decay rate with respect to the first quark energy. The yellow dash 
line represents the virtual contribution from three body decay which is rescaled by a 
factor of 3.24 x 10 2 , while the blue line represents two body decay result. In both 
plots, the initial KK gluino mass is Pg = 1.15 TeV and we use g§ = 0.108 TeV -1 / 2 . In 
the left one, we choose z uv = 10 -3 TeV" 1 , z m = 80 TeV" 1 , c = 0.5, fi g = 0.4 TeV 
and fj,f — 0.5 TeV; In the right one, we choose zuv — 10 -3 TeV -1 , zm = 100 TeV -1 , 
c = 0.5, fig = 0.3 TeV and fi f = 0.4 TeV. 

masses of the initial and final gluinos: 

E <min = 0, (43) 
1 / (p~) 2 \ 

The contribution to the differential decay rate as a function of the outgoing quark 
energy from the virtual 3-body process can be compared with the real 2-body contri- 
bution by summing over all final KK gluinos and outgoing anti-quark energies which 
are kinematically allowed. The result is shown in Fig. |5j The virtual 3-body process 
is the cut of a two-loop diagram while the real 2-body process is the cut of a one-loop 
diagram, so we expect that the 3-body decay rate should be suppressed with respect to 
the 2-body decay rate by a 4D loop factor as long as the theory remains perturbative. 
The shapes of the two contributions are also similar, but with the 3-body contribution 
being slightly harder. Thus we can conclude that in this case, it is reasonable to use 
the narrow-width approximation to calculate the energy distributions of the visible 
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P,z IR dr/dEu° L P 9 z IR dr/dEu !, 




Eu° L (Tev) 



0.1 0.2 0.3 0.4 0.5 . 1 0.2 0.3 0.4 . 5 .6 

Figure 6: 2-body gluino differential decay rate as a function of the outgoing quark 
energy E u o in TeV. We have taken in this example Zuv — 10~ 3 TeV -1 , Zjr = 80 TeV -1 , 
c = 1/2 and g§ = 0.108 TeV -1 / 2 . In the figure on the left, we fixed the initial gluino 
KK-mass at p g = 1.26 TeV and its mass-gap to /x s = 0.3 TeV, and vary the squark 
mass gap by /// = 0.40 , 0.43 , 0.46 , 0.50 , 0.56 TeV. As can be seen from the figure, the 
peak position decreases in magnitude and shifts towards larger values of E u o . In the 
figure on the right, we fixed \x g = 0.4 TeV and \x$ = 0.5 TeV, and vary the initial gluino 
KK-mass by p g = 0.83 , 0.97 ,1.15 , 1.29 , 1.52 TeV. In this case the peak increases in 
magnitude with increasing p and its position remains roughly constant as a function 
of E..o . 

a L 

particles coming from the continuum decays by treating each decay step going to real 
states that subsequently decay independent of the details of the previous decay. 



4 Phenomenology of the Continuum Superpartner 
Decays 

The comparison of the differential decay rate calculations for a real 2-body process 
and a virtual 3-body process in the previous section, shows that it is reasonable to 
calculate the energy distributions of the visible particles coming from the continuum 
superpartner decays using the narrow-width approximation. This greatly simplifies the 
phenomenological study of continuum superpartners. As mentioned earlier, the contin- 
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Figure 7: Same as Fig. [6] for different sets of energies and mass gaps. In the plot 
on the left, we fixed p g = 2.4 TeV, and \i g = 0.2 TeV. We vary the squark mass 
gap by jif = 0.36 ,0.40 ,0.43 ,0.46 ,0.5 TeV. As can be seen, the peak decreases in 
magnitude and slightly moves towards larger values of E u <^. In the figure on the right, 
we fixed \x g = 0.2 TeV and fif = 0.3 TeV, and vary the initial gluino KK-mass by 
p g = 1.24 , 1.62 ,2.01 ,2.21 ,2.4 TeV. In this case, the peak magnitude increases and 
its position remains roughly constant with respect to E u o . 

uum states tend to decay to other continuum states that are nearby in invariant mass, 
and so the ordinary particles tend to be emitted with soft energies. This is the behavior 
that is expected in a CFT [814TT]. Thus if the decay chain starts fairly high up in the 
continuum, then there is usually a long decay chaining with an approximately spherical 
distribution of energy [TTJ. In this section, we examine how the energy distributions of 
the visible particles depend on the parameters of the theory and the process. 

As the energy of the emitted quark E u o increases (or what is the same, the squark 
mass q n decreases), the phase space of the final particles increases while the vertex 
c(p n , Qm) decreases. It is the competition between these two factors that makes the 
differential decay rate have its maximum at an quark energy around E u o^ max ~ (a*/ - A*g) 
as can be seen in Figs. M and [7J As a consequence, for a high-energy initial state 
(p g ^> [if, fig), we expect order p g /(fif — p, g ) particles coming out a decay chain. For 
fixed #4, c, zir, and zjjy, there are only three parameters: the initial gluino energy 
Pg and the two mass gaps ///, fi g . Since there should be little dependence on z IR and 
ZuV) as l° n g as they are far away from the energy scale interested, physical quantities 
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Figure 8: u L>™ ax vg> j±g anc j (j ) _ ZlR ^_^_\^ max vs _ N—Bl evaluated at the peak 
positions. We use z uv = 10~ 3 TeV" 1 , z m = 80 TeV~\ c = 0.5 and g 5 = 0.108 TeV~ 1/2 
and fix p g = 2.40 TeV, \i g = 0.2 TeV, while the squark mass gap is varied from 
fMf — 0.3 TeV to /if — 0.75 TeV. The plots display the fitted functions. 

can be expressed as functions of the two dimensionless ratios of p g , {if, and /i 9 up 
to an overall normalization. In Fig. |H] we show the dependence of the peak position 
and magnitude of the differential decay rate on the initial energy p g and the mass gap 
difference fif — fi g . From these figures we can determine an approximate behavior for 
the peak position, i? u o max in the differential decay rate as a function of fif — fi g and 
p g . We find that it can be parameterized by an exponential decaying function with an 
overall normalization A and a numerical exponent B: 

4»,max = MVf ~ P>g) e P ' 9 ■ ( 45 ) 

The numerical values of A and B are obtained by fitting numerical data points, as shown 
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in Fig. |H1 where we find that A ~ 1 and B 3/2. Using the plane wave approximation 
for the wave functions of the incoming gluino and outgoing squark, the coupling squared 
is proportional to \c(p m ,q n )\ 2 oc z $ jjtp^ ) 2 ■ Then one can substitute Eq. (l4"5l) into 
Eq. fl24l) . using AE u o l « k/zir, to obtain an approximate analytical expression for 
(dT/dE u o )\e o as a function of /i/ and fi g : 



ii~ ,max 

dT 

PgZlR 



dE..o 



-2B 



(M/-Mg) 



C-gi 



e 



2 (log(l/(2/x g z^))-7) 
4vr 2 (1(^(1/(2,//^)) - 7 ) 



(46) 



the overall normalization C compensates for the inaccuracy of the plane wave approx- 
imation, and we have written everything in terms of 4D quantities. The dependence 
on zuv is only logarithmic and therefore mild. To give an idea of the validity of the 
approximation, we plot against numerical data in the second plot of Fig. [HJ We see 
that ( 1451) and ( 1461) provide reasonably good estimates for the functional dependence. 



5 Conclusions 

Supersymmetry is one of the best motivated scenarios for new physics beyond the 
standard model at the TeV scale. For the past two decades it has been intensively 
searched for. Currently, the experiments at the LHC have placed very strong limits on 
the masses of the squarks and gluino to be above ~ 1 TeV in the standard scenario [TBT 
1T8] . This means that either the superpartner spectrum is unnaturally heavy or the 
superpartners decay in unusual ways which escape the standard SUSY searches. As we 
showed in this paper, if the superpartners have continuous spectra, they tend to have 
long decay chains and produce many soft SM particles. This is a challenging scenario 
at the LHC because the soft particles may not pass the experimental cuts, and the 
signals could be buried in the QCD backgrounds. It may require a more specialized 
study to search for this kind of signal. 

Naively the phenomenological studies of continuum superpartners at colliders may 
appear to be formidable as the usual calculation techniques have only been developed 
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for particles. Here we have seen that these methods can still be used in continuum 
calculations. In particular, we showed that the narrow-width approximation is still 
generally valid in the perturbative decay processes that we are interested in. This 
greatly simplifies the calculations because processes with long sequences of decays can 
be divided into individual steps involving real states and each step can be calculated 
independently. The easiest way to perform the calculations is to introduce a regular- 
izing IR brane in the 5D picture which transforms the continuum into discrete KK 
modes, then one can carry out the calculations as in the usual particle case. As long 
as the KK modes have reasonably fine spacings, the results are basically independent 
of the position of the IR brane. In this way, continuum superpartners can also be 
implemented in the usual collider simulation tools in order for more detailed studies to 
develop new strategies to search for such kinds of exotic collider signals. 
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A Equivalence of the Continuum and KK Limit 

It is interesting to compare our expression for the propagator. Eq. f )39|) . with the more 
familiar KK-representation of the 5D propagator which can be expressed as sum over 
the different KK-level propagators: 



where fL(m n ,z) is the normalized left-handed squark wavefunction. This propagator 
corresponds to a real particle when q 2 = vr? n (aka going "on-shell") for each particular 




(47) 
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KK-mode n. For momenta q close to m n , we can use the approximation, 



q 2 — m 2 + ze ' " q 2 — m 2 

where V denotes the Cauchy principal value. As the mass difference between adjacent 
KK-levels n and m = n — 1 tends to zero (the continuum limit), Am n = m n — m m — > 0, 
the sum over KK-levels in Eq. ([471) becomes an integral in the complex plane over 
"tightly squeezed" resonances which form a branch cut in the limit zir — » oo. Thus, we 
can identify the real ("on-shell") contribution of the continuum propagator, Eq. (f3"9"j) . 
by matching the continuum limit of Eq. (j4"Tj) with Eq. (1391) . For that purpose, we 
consider Eq. (jUJ in the limit /i/ < g < 1/zjjv- Let us take z — z' — zjjv and using 
Eq. f[T3"j) and Eq. ffTTI) . in the case c = 1/2 for simplicity, we find that: 



h(m n ,z uv ) 



2 



— >■ 



. M(0, 0, 2?m n ^) 

- u/ / n n o- rW(0, 1, 2«m^ C /v / 

VK(0, 0, 2im n z C ry) 



^ i (49) 

*™ v/^av ((log (a.) + log(^ y ) + 7 e) 2 + f ' 

By identifying Am n = n/z IR and using Eq. (]4"5]). we find that: 

ImfP^g, ^[/y)] -> /i(m nj ^y) 7r5(g 2 - m 2 n ) 



- ^ -o — ^ = M— ]■ (•><» 

4 



g 2 ^y ((log (I) + log{z UV ) + lE f + ^ 



One can also see the equivalence between the KK picture in the limit zir — >■ +oo and 
the continuum picture by comparing the differential decay rate for large zjr in the KK 
picture versus the continuum case where one still keeps the particle behavior for the 
initial gluino state, but replaces the final state squark phase space, by an unparticle 
phase space [2]. 

For that purpose, let us calculate the continuum squark phase space for c 7^ 1/2 and 
c = 1/2. In order to get the phase space of the squark final state, we calculate first the 
spectral function p(p 2 ) = 2 • lm(iA(p 2 )), where A(p 2 ) is the squark correlator whose 
expression was calculated in Ref. pp. As was done in previous calculations, we are 
interested in the case when the momenta involved are much bigger than the squark's 
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mass gap, p 3> Using the expressions of the correlator in the limit when pzuv 1 
we then find that for -1/2 < c < 1/2, 

2 _ 2B{p 2 ) 

where 

.2 e i- 2c 2 - 1+ V(P 2 - /4)~ 1/2+c cos(tt(-| + c))r(l - 2c)r(c) 



A(p : 



P 



i - 2c r(2c)r(i - c) 



2-^gV - tf )-^sm(n(-± + c))r(l - 2c)r(c) 

" r(2 C )r(i - C ) ' {bI) 

and that for c = 1/2, 

P(P 2 ) = — • (53) 

( 7s + log(f v /p 2 -^ ))2p2 + ^ 2 

The continuum squark phase space is then given by, 

d$ g = % )% 2 -/4)p(p 2 ), (54) 

where 8(p) is the Heaviside step function. On the other hand, the phase space for the 
outgoing quark is just the usual particle one, 

rf$ e = 2n6(p°)5(p 2 ). (55) 

With all this information, we are ready to calculate the 2-body phase space integral, l 
Let us name the momentum of the incoming gluino as p gi and of the outgoing quark 
and squark momenta as p u and p u . Then the 2-body phase space integral is given by 

n 2 = f ^)^\p- g -p,-p u )2^{pl)6{p^ (56) 

We can perform the integral over the part of the quark phase space using that, 

(57) 



o xr2\arO\dPu f d Pu 



(2vr) 4 J (2vr)32£ u 



P°u>0 



4 Even though the vertex depends on the momenta of the particles, it only ends up depending on 
the quark energy E u and the initial gluino energy . 
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Figure 9: We use z uv = KT 3 TeV -1 , z IR = 80 TeV -1 , c = 0.5 and g 5 = 0.108TeV~ 1/2 
and fix p g = 2.40 TeV, fi g = 0.2 TeV, and /i/ = 0.3 TeV. The plots compare the 
continuum method (red dashed curve) with the KK mode method (blue curve), they 
are almost identical. 

Using the 4-dimensional delta function to perform the integrals over d 4 pn and in the 
center of mass frame (CM) of the gluino, we can trivially perform the angular integra- 
tion, d 3 p u = EldE u dQ, so that we obtain, 

n 2 = ^^- u P((pD 2 - 2P° S E u )e((pl) 2 - 2p\E u - fi)dE u . (58) 

Once we have the expression for the phase space integral, we can write the differential 
decay rate as, 

\M\ 2 

dv = L_^n 2 . (59) 

Averaging over the initial and summing over the final spin and color, we find that 

\M\ 2 = 2\d(p § ,E u ,c)\y § E u ti[t a t a }/8 , (60) 

where 

n 1 c(p~, E u , c) 

zjy /(?> z uv) 

is the 4D effective vertex obtained from integrating over the profiles in 5D, divided by 
the normalized squark profile evaluated at the UV brane, and multiplied by to be 

z uv 

dimensionless, with c{p° g) E ul c) shown in Eq. (|22|) . Thus we arrive at the formula, 

^ ~ El |^,K,c)| 2 p((p°) 2 -2^)e((^) 2 -2^- / ,J). (62) 



dE u 4tt 2 
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In Fig [9] we plot the KK differential decay rate, Eq. f l24j) . and the continuum decay 
rate, Eq. ( 162]) . as a function of the final quark energy E u . They are in excellent 
agreement for large enough zir. 
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